Proposal of a full Bell state analyzer for spin qubits in a double quantum dot 
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We propose a scheme for full Bell state measurement of spin qubits in a double quantum dot. Our 
scheme consists of Pauli spin blockade measurements and biaxial electron spin resonance. In order 
to eliminate the average of the Zeeman fields, the double quantum dot is designed so that the Lande 
g-factors of first and second dots satisfy gi = — g2 with the use of g-factor engineering. Thus, we 
can swap one of the three spin-triplet states for the spin-singlet state without disturbing the other 
states. Our study shows that the sequential spin-to-charge conversions enable us to implement the 
full Bell state measurement of electron-spin qubits. 
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\ Bell state measurement is a key element of quantum 
information science and technology. It generates an en- 
tanglement state of qubits, and plays an important role 
in quantum teleportation^, entanglement swapping^ and 
quantum key distribution^. On the other hand, much ef- 
fort has been devoted to realize solid-state qubits, which 
are promising candidates for quantum memories and 

^ quantum gates. Bell state analyzers for these qubits are 
then highly desired to explore the potentiality of solid- 

\ based quantum devices. 

I In 2005, Engel and Loss proposed a partial Bell state 
analyzer for electron spin qubits that determines the par- 
ity of two qubits^.. However the partial analyzer re- 
quires some initial source of entanglement for the full 
Bell state measuremenl4. In this Rapid Communication, 
we propose a protocol and physical implementation for 
full Bell state measurements of spin qubits in coupled 
quantum dots (QDs). Our scheme consists of Pauli spin- 
blockade measurements^"— and biaxial electron-spin res- 
onance (ESR) ''"^'^. We show that the sequential spin-to- 
charge conversions and total spin rotations enable us to 
implement the full Bell state measurement if the Lande 
g-factors of the first and the second QDs are designed to 
satisfy gi = — g2 with the use of g-factor engineering. 

The system we have in mind is vertically or later- 
ally coupled QDs as shown in Figs. [Ija) and [Hb). 
We employ the so-called Hund-Mulliken model where 
the orthonormal state is defined as $1(2) (r) = (^1(2) ~ 
1 5V2(i))/\/l ~ SFg -f g'^^^. Here F is the overlap inte- 
gral of the orbitals in first and second QDs 951,2, and 
g = {1 — \/l — F^) /F. The relevant two-electron states 
arc the (1,1) charge states, {jS"), |To), |r+), |r_)}, and 
the spin-singlet (0,2) charge state, | (0,2)5). The la- 
bel (m, n) refers to the electron number confined in the 
first and the second QDs, and the Bell basis set is com- 
posed of the four (1,1) states, i.e., one spin-singlet state, 
\S) = [| tl> - Ut>]/V2, and three spin-triplet states, 
iTi)) = [|n) + i;t)]/v^ and |r±) = [|tt) ± IU)]/v^. 
Throughout the work, we neglect the spin-triplet (0, 2) 
state because it is energetically inaccessible^. The Hamil- 



tonian of the double QD is 

^ ~e|(0, 2)5) ((0, 2)S\ + V2T [|(0, 2)5) (5| + H.c] .(1) 

Here e and T represent the gate-controlled detuning en- 
ergy and the hopping integral both of which are functions 
of F, and H.c. stands for Hermitian conjugate. 

When a magnetic field B is applied to the double QD, 
the Zeeman energy is described by the Hamiltonian //z — 
^^^-^ 2 hj • Sj/fi., where Sj is the spin-1/2 operator for 
the j-th QD and = — gj^eB with /^b being the Bohr 
magneton. Here the Lande g-factor is defined by gj = 
(<I'j|g(r)|<I>j) since we treat spatially varying g-factor be- 
low. Let us introduce the two-spin operators s = §1 + §2 
and (5s = §1 — §2, and the average (inhomogeneity) of the 
Zeeman fields h = (hi + h2) /2 (5h = (hi - h2) /2). In 
terms of the Bell basis, (|5), ITq), 1?+), |r_)), the Hamil- 




FIG. 1: (color online) Schematic view of (a) vertically and 
(b) laterally coupled double QDs. Both the QDs are located 
between the heterostructure barriers, and are then focused to 
small areas with the use of the surface metallic gates. We as- 
sume that the g-factors in the two QDs are designed to satisfy 
the condition gi = — g2. For the vertical double QD, such a 
condition can be achieved, e.g., by adjusting the thickness and 
the alloy composition of each quantum well. For the lateral 
QDs, we can achieve the condition by controlling the equilib- 
rium electron positions in individual QDs when the top and 
bottom barrier layers have different g-factors. We take the 
z-axis as the direction of the static magnetic field, and apply 
oscillating magnetic fields in the x- and y-directions common 
to the two QDs. The ac fields flip the electron spins confined 
in the QDs when the ESR condition is satisfied. 
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tonian of the Zeeman energy is expressed as 
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When the average of the Zeeman field is zero, i.e., h = 0, 
one can see that Eq.([2|) reduces to 

Hz = Sh'^\S){To\ - i6hy\S){T+\ ^ Sh^\S)(T- \ +U.C.. (3) 

Thus we can sequentially swap one of the triplet states 
for the \S) state without disturbing the other states by 
controlling the inhomogeneity 6h.. 

We assume that the system is designed so that the 
condition gi = — g2 is satisfied using g-f actor engineer- 
ing to eliminate the average Zeeman field. In general, an 
electron g-factor is a function of band gap and spin-orbit 
interaction in a semiconductor^^,. Therefore, it is possi- 
ble to adjust the electron g-factors by changing the alloy 
compositionsi^iii or well thickness to nearly zeroi^iia. In 
addition, the structurally-defined g-factor is affected also 
by the penetration of the wave function into the bar- 
rier layers^. Then it is further controlled by shifting 
the equilibrium electron position between the two layers 
with different g-factors. Indeed, the g-factor control2£ 
and fine tuning across zero™ have been achieved using 
a heterostructure or a quantum well. The g-factor en- 
gineering is applicable also to single electron in a QD. 
We have actually fabricated the g-factor engineered QDs, 
and electrically confirmed that the g-factor is as small as 
designed (g = 0.05)22.. 

Transition between up and down spins of a local- 
ized electron is accomplished using ESR2^. Here we 
take the z-axis in the direction of the static magnetic 
field, i.e., Bq = (0,0, i?o)- Thus the ac magnetic field 
Bac(i) = Bicoswf in the x-y plane makes the electron 
spin fiip when it is resonant with the spin precession. 
However, it may be challenging to produce strong and 
oscillating fields in more than one direction. There is an- 
other experimental technique, called electric-dipole spin 
resonance, that uses the electric gate^rJ^. The effective 
ac field is generated by the cooperation of the oscillating 
electric field Eac(i) = Eicoswt and spin-orbit interac- 
tion in the two-dimensional (2D) QDs^iiiiH. Within the 
linear momentum regime, the effective field becomes^iii^ 



Bcff(t) = Bq 



A. 



(4) 



where hwd (Id) is the confinement energy (lateral size) of 
the 2D QDs, and x' = {x + y)/V2, y' = {x-y)/V2. The 
spin-orbit length is defined as X± — h/me{l3 ± a), where 
nic is the electron mass, and a {(3) is Rashba (Dressel- 
haus) spin-orbit coupling constant^"*. 

We show below that manipulating e and Bac(i) en- 
ables us to conduct the full Bell state measurement. Ini- 
tially the detuning energy is set to be large negative, e = 
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FIG. 2: (color online) The sequence of the Bell state mea- 
surement we propose is shown schematically. The solid lines 
indicate the energy costs to add an extra electron to the first 
or the second QD. The dotted line with a filled circle indicates 
the energy levels of the (0,1) charge state, (i) Initially the de- 
tuning energy is set to be a large negative, e = eo <C —y/2T. 
(ii) Applying gate voltage to the first QD (depicted in left), 
we apply the spin blockade regime (e = ei > V2T). The 
electron in the first QD can tunnel to reservoir when the two- 
electron state is in \S). (iii) If no electron tunneling is de- 
tected, the system is brought back to the Coulomb blockade 
regime e = eq. We apply the effective ac field B^^{t) for the 
time T-ESR on both the QDs, and then we wait for the time 
to- This swaps |r_) for \S}. (iv) We apply the spin blockade 
regime to the QDs, and determine whether the two-electron 
state was initially in |r_). If the electron tunneling is not 
detected again, (v) we flip |r+) into l^) by B^^{t), and (vi) 
determine that the system was in |r+) or \To). 



eq <C —V^T, so that the (0, 2) charge state lies far above 
the (1,1) states as shown in Fig. [SJi). Then, any two- 
electron state can be expressed by superposing the four 
Bell states, i.e., \^p) = c,|5) + co|To) + c+|r+) + c_|T_). 
For Bell state measurement, it is necessary to perform 
projection measurements of \tp) to all the Bell states. 

We begin with the projection measurement to \S). We 
adjust gate voltages on the QDs to apply the spin block- 
ade regime, i.e., e = ei > \/2T as in Fig. [2][ii). When 
the gate voltage is swept adiabatically for the interdot 
tunneling \/2T, the singlet component of IV') segues into 
I (0, 2)S) through an avoided crossing. Note that the gate 
sweep should be nonadiabatic for the inhomogeneity of 
the Zeeman field SHq = — gi/iB^o so that the mixing of 
\S) and I To) can be inhibited2^'2S. Setting the Fermi level 
of the drain to the energy of |(0, 2)5'), the interdot tun- 
neling is detected as a current. Because the triplet (0, 2) 
state is neglected, the detection of the electron tunneling 
determines that \S) is postselected^"— . 

In the null result case, the postmeasurement state be- 
comes IV') = N[co\To) + c+|T+) + c_|T_)] with normal- 
ization factor N = 1/(1 — |csp)^^'^. We thus move on to 
the next stage, i.e., projection measurement to |T_). To 
measure |T_) via the Pauli spin-blockade measurement 
described above, we should swap |T_) for IS"). Therefore, 
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we apply the oscillating magnetic field B^^{t) common to 
both the QDs after bringing the energy detuning back to 
e = eo as shown in Fig. [2] (iii) . In the rotating frame of 
reference, the two-electron spin state is expressed as 



my) = exp 



1 



j2igjf^BBos^M\m)- (5) 



When the resonance condition, huj = \gj\fj,BBo, is satis- 
fied, the electron spins rotate in the rotating frame as 



If), ^ n^mr), = cos ^ir>, + ism^\v),, 

-> 7^|(^)|r>, = cosli;-^), + zsin||r)„ 



(6) 



where 9j = Qjt is the rotation angle with Rabi frequency 
ilj ~ gjfisBf /2h. Note that the spins in the two QDs 
rotate in opposite directions. If the burst time for the 
ESR tesr satisfies the condition that 9i{tesr) = 7r/2, 
|Ti) is transferred into {S") without disturbing the other 
two states. In the laboratory frame of reference, however, 
the coherent oscillation between 15) and \To) evolves due 
to Then, after the ESR sequence, we should let the 
oscillation evolve for the time to which is determined by 
the condition that 5/io(io + ''"ESRj/fi- = (mod 2tt). As 
a result, the two-electron state is transferred as 1-0) 
7^^(^o)7^"(TESR)|V') = N[co\To) + c+\T+) - ic-\S)]. We 
then perform the Pauli spin blockade measurement to 
determine whether |T_) is postselected (see Fig. [^Kiv)). 
When the null result is gained again, we swap |r_|_) for 
IS") by By^{t), and determine which of \T+) and |To) is 
postselected, as shown in Figs.[2](v) and[2jvi). 

In the preceding discussion, we have assumed that 
h = 0. However, in the real system, there should be 
a small mismatch of the absolute values of the g-factors, 
i.e., 5g = — (gi+g2) 7^ 0. This produces a small but fi- 
nite h in Eq. which will reduce the fidelity of the 
swapping operations. When the ESR frequency is set 
to w = \ga\HBBo/h with go = (gi - g2)/2, the mis- 
match of the ESR frequency in each QD is given by 
6ujj = —{—iy6giiBBo/2h. Here we assume that the fre- 
quency mismatch Sui and the Rabi frequency in each QD 
rij are much smaller than \gi\fiBBQ, and thus the rotating 
wave approximation is justified. The ESR Hamiltonian 
of the double QD becomes^^ 



ESR 



E 
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(7) 



in the rotating frame. For instance, let us estimate the 
fidelity of the swapping between |T_) and \S) states. We 
apply the oscillating field B^^{t) to both the QDs for 
''^SR = 7i'/4r2o, where J7o — ^bB/AH is the average of 
the two Rabi frequencies with B — yj ((5gi?o)^ + (goSf )^. 
The resultant fidelity is calculated as 



l(5■■|7^^(rESR)|T:; 



IcosOK 1 + C 



{bgBlf 



B^ 



where tanO = SgBo/goBf and SgBf B. To achieve 
the fidelity > 99%, we have to control the g-factors so 
that Sg/go < 0.01 when the magnetic fields {Bo,Bi} = 
{41mT, l.QmT} in Ref. are employed. 

In addition to the g-factor mismatch, there are a vari- 
ety of sources for measurement errors. For instance, the 
spin-orbit interaction causes spin fiip through electron- 
phonon interaction a ^^ ' ^'' and total spin rotation in inter- 
dot hoppin g^^i^^ . However, the most relevant difficulty 
is hyperfine interaction with host nuclei. The hyper- 
fine Hamiltonian is given by Hbf = '^j=i,2^nj ■ Sj/h. 
The root mean square of the field reaches |(hnj)i.i„s| 
10~''meV in a QD containing unpolarized N = 10^ nu- 
clear spins^i^i^. Since the hyperfine field cannot be 
controlled by the electron g-factors, the fluctuating nu- 
clear spins shorten the spin dephasing time T2 and de- 
crease the fidelity of the singlet-triplet swappings. 

It is possible to prolong Tj* by preparing the nuclei in 
an approximate eigenstate of Hbf before we conduct the 
Bell state measurement. Here we consider the method 
for dynamic nuclear polarization^^''^^ in the QDs with 
the g-factors gi = — g2. Note that it is desired to pre- 
pare the nuclei so that hn ~ because the large aver- 
age Zeeman field is inimical to our Bell state analyzer. 
The energy diagram in Fig. [3Ka) shows that there exists 
a state crossing E — 0^— which can be expressed as 
|0) = sm(j)\To) - cos 015(0, 2)) with tan0 = Sh^y^T 
at zero detuning. The hyperfine field becomes relevant 
around E = 0. The component parallel to the exter- 
nal field, hi, lifts the degeneracy between |tt) and Hi), 
whereas the perpendicular part opens the gaps nearby 
the crossing point (see Fig. [3Kb)). We first prepare the 
system in |(0,2)S') with detuning e = eo, i.e., position 
A in Fig. [3Ka). By sweeping the detuning to ei adia- 
batically and bringing it back to eo nonadiabatically, we 
apply the two-electron state in point C. Then, by driving 
the system adiabatically along |0), the flip-flop process 
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FIG. 3: (color online) (a) Energy diagram for the relevant 
states. From initial point A with S{0,2), we can reach point 
C by sweeping the detuning e adiabatically (A— >-B) and then 
nonadiabatically (B— >C). (b) Enlarged view of the region in- 
dicated by square in panel (a). The hyperfine field allows for 
flip-flop process (\To) — >■ |tt)) when ftil > 0. Then driving the 
detuning up adiabatically from C can reduce the total hyper- 
fine field /ill. For < 0, flip-flop process (|To) — >■ lii)) can 
reduce \hn\ in a similar way. 
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between |To) and |tt) can occur when /iJl > 0. This pro- 
cess changes the average hyperfine field along Bq, i.e., 
reduces /lil. When ht < 0, on the other hand, it is |||) 
that |0) surges into. As a result, cyclic repetition of these 
flip-flop transitions leads to a hyperfine field with small 
h„ and large Shn- Although this method seems complex, 
any of the processes has been already achieved in exper- 
iments using a double QD with gi = g2^i^. 

As for the time needed in experiment, the burst time 
is 27 ns to rotate a single spin by 7r/2 for a QD with |g| = 
0.4 and Bi = 1.9 mT^. In each spin-to-charge conversion 
it takes r ^ 1 ns to sweep e for ^ 1 meV, adiabatically 
for v^T ~ 10-2 meV^>26 but nonadiabatically for - 
10"'^ meV— . Such a condition enables us to inhibit the 
S-Tq rotation due to |(5/i§| and the g- factor modulation 
by the gate control^^. The usage of vertical QDs with 
large coupling, e.g., ^/2T ~ 0.3 meV§ can decrease r 
and then freeze the rotation further. Besides, a resonant 
tunneling to the drain is estimated to take 1.6 ns in Ref.i. 
Therefore, if the dephasing time Tj* reaches ~ 1 /xs as in 



the double QD with gi = g2^, one can verify our scheme 
in experiment. 

In conclusion, we have proposed a full Bell state ana- 
lyzer for spin qubits in a double QD that are designed to 
have g-factors with opposite signs, i.e., gi = — g2. We ex- 
pect that the proposed analyzer enables various devices, 
for instance, a quantum repeater that is essential for 
ultra-long quantum communication^. Indeed quantum 
state transfers from 'message carrying' photon qubits to 
"storage" matter qubits, e.g., atom ensembles^i^^, ^^P 
nuclear spins^^ and electron spins in a semiconductor— 
were demonstrated. Although the transfer to electron 
spins were confirmed in ensemble measurements in a 
quantum well^* , the mechanism works for single electron 
in a g-factor engineered QD with slight modification^—. 
Combined with these results, our scheme may contribute 
to constructing quantum information networks as well as 
processor units. 
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